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Abstract 

In this communication, the approach of phenomenological universalities of growth 
are considered to describe the behaviour of a system showing oscillatory growth. 
Two phenomenological classes are proposed to consider the behaviour of a sys¬ 
tem in which oscillation of a property may be observed. One of them is showing 
oscillatory nature with constant amplitude and the other represents oscillatory 
nature with a change in amplitude. The term responsible for damping in the 
proposed class is also been identified. The variations in the nature of oscilla¬ 
tion with dependent parameters are studied in detail. In this connection, the 
variation of a specific growth rate is also been considered. The significance 
of presence and absence of each term involved in phenomenological description 
are also taken into consideration in the present communication. These proposed 
classes might be useful for the experimentalists to extract characteristic features 
from the dataset and to develop a suitable model consistent with their data set. 
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1 Introduction 


In the field of science and technology, it is found that many characteristic pat¬ 
terns have been discovered that are remarkably similar in nature, though they 
are related to totally different phenomenologies [31S H [S]. It is because of 
the fact that the background mathematical formalism remains the same. There¬ 
fore, an universality in terms of mathematical formalism is observed, in which 
a mathematical equation is showing universality in the sense that the variables 
involved in the equation may indicate a scientific discipline. The opposite may 
also be true in the sense that some related effect might be found in nature for 
every mathematical niche, that is generally observed in biology. However, this 
might be found in the other field of application. For example, any symmetry 
class, predicted mathematically, is found to be available in nature. Therefore, a 
search for universality in different phenomenologies is very important from the 
point of view of applications. A cross-fertilization among different fields might 
be accelerated when different fields share a common mathematical formalism or 
concept. One of such approaches is the classification scheme of phenomenologi¬ 
cal universalities [5]. 

In the phenomenological approach, the universality class at the level n = 0 
corresponds to Malthusian or autocatalytic processes. At the level n = 1, it 
produces Gompertz law of growth [7] that was being used for more than a cen¬ 
tury to describe growth phenomena in diversified fields. A variation of the class, 
termed as involuted Gompertz function, is used to describe the evolution and 
involution of organs undergoing atrophy [S]. The class corresponding to the 
level n = 2 leads to the growth equation, proposed by West et. al. to describe 
different biological growth patterns [iiiniiiiiiiiiis]- The classes correspond¬ 
ing to n = 1 and n = 2 have been applied in diversified and unrelated fields 
[n [13 [m [H]. An extension of this approach, termed as complex universal¬ 
ity, was characterized [TH] and applied to explain the concurrent growth of two 
phenotypic features |19j . A vector generalization of the phenomenological uni¬ 
versality is used to describe the interactive growth of two or more organisms in 
a given environment [201121] . 

The formalism of phenomenological approach [6] of classes n= 0, 1, 2 leads to 
monotonous growth of the system with no scope of oscillations, though they are 
an essential feature of different complex systems of nature. In some cases they 
are not observed, might be due to the fact that they were overlooked, or masked 
by large experimental errors or missed due to the inadequate range of variation of 
the variables being measured. Apart from these factors, oscillations are observed 
in different types of complex systems [131131 [Ml 1131 HSl HZl UHl IMl 13011311133 • 
The existence of volume oscillation, of relatively large amplitude, in the growth 
of multicellular tumour spheroids was experimentally found [331134] . A varia¬ 
tion of the phenomenological approach in the complex field, instead of real field, 
termed as complex universalities is applied to analyze two correlated features of 
the growing oscillatory systems [18] . In the same framework, concurrent growth 
of phenotypic features are also explained and the correlation between pheno¬ 
typic features is derived by Barberis et. al. [H! .One of essential features of the 
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complex universalities approach is that the growth or decay in the amplitude is 
unavoidable [18]. Therefore, it is not possible to explain the oscillating behaviour 
of a system with constant amplitude. 

In the phenomenological approach, it is expected that the values of the coeffi¬ 
cients are governed by the particular growth mechanism involved in the system 
|5]. The biological and environmental constraints might induce a new growth 
processes or cease a growth mechanism already existing in the system. There¬ 
fore, the inclusion or withdrawal of a growth mechanism from a system may be 
considered as deviation from the normal growth process due to some constraints. 
The allometry based biological growth is represented by the phenomenological 
class 172 [^. Normal growth in a biological system shows a monotonous growth 
process by nature and attains a saturation level in course of time. A deviation 
from normal growth in terms of oscillation may be treated as adaptation or 
withdrawal of a new growth mechanism represented by different phenomenolog¬ 
ical terms other than a + ba'^ [B] . The growth of a tumour may be treated as 
a deviation from the normal growth due to some biological and perhaps envi¬ 
ronmental constraints. The coexistence of tumour growth along with normal 
growth of tissues within the same biological system and the beginning of tumour 
growth from normal growth may be treated as an indication of such deviation. 
The growth of a tumour shows such oscillatory nature in terms of its volume 
with relatively large amplitude |331 134). This communication is motivated by 
these phenomenological observation and aimed to describe phenomenologically 
oscillatory behavior of a system. 

In the present communication, phenomenological approach based on the formal¬ 
ism presented by Castorina et. al. [6] is considered to describe the behaviour 
of a system showing oscillation. Here, we propose two different phenomenolog¬ 
ical classes. One of them shows oscillation with constant amplitude (class-I). 
The other shows growth or decay in the amplitude of oscillation (class-II). The 
deviation from the proposed phenomenological class-I in terms of addition of a 
phenomenological term with the specified condition, might lead to a damped 
oscillatory growth (class-II). The emphasis is given to find out the role of coeffi¬ 
cients involved in the phenomenological description of a dynamical system. The 
nature of variation of specific growth rate with different parameters is also stud¬ 
ied. Absence of each phenomenological terms in the proposed phenomenological 
classes is also examined. Even a linear growth feature of the system may also be 
observed in absence of a phenomenological terms related to normal growth pro¬ 
cesses. In each case, the role of the phenomenological coefficients are examined 
in detail. The phenomenological term responsible for damping in the proposed 
class is also identified. This analysis is focused to help the experimentalists to 
decide whether their dataset, showing oscillatory nature, might be considerably 
expected to fit with these phenomenological classes. 

The paper is organized as follows: In Sec. II, we would propose two phenomeno¬ 
logical classes presenting oscillatory nature. In this connection, the classification 
scheme of phenomenological universalities proposed by Castorinal et al. [6] is 
discussed in brief. Different aspects of proposed phenomenological classes would 
be considered in Sec. III. Finally, we would conclude with our results in Sec IV. 
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2 Phenomenological approach to oscillatory growth 


Following the formalism of Castorina et al. [6], a string of data show¬ 

ing a particular phenomenon that can be described by an ordinary differential 
equation, is represented in phenomenological approach in a following manner [5], 

^ = a{t)Y{t) (1) 

where, a{t) stands for specific growth rate of the variable Y{t). Equation ([T]) 
can also be expressed in terms of adimensional variable as. 


= a{T)y{T) ( 2 ) 

Where, r = a{0)t, y{T) = Y{t)/Y{Q) and a(T) = a(t)/a(0) are adimensional 
time, adimensional variable and adimensional specific growth rate respectively. 
It is assumed that a(0) = 2/(0) = 1. An adimensional variable (z) is introduced 
in this connection so that, 

z = lny (3) 

The time variation of z is described by the following relations, 



(4) 


and. 


ip{a) 


(Pz 


Again, ip{a) can be expressed in terms of power series as. 


(5) 


( 6 ) 

0 

From equation ([2]), ([3]) and ([5]), the following expression can be derived. 


z = 


adr -I- constant 


(7) 


Different forms of would generate different types of growth equation. Each 
term in the right hand side of the equation ([5]) represents different types of 
growth mechanism involved in the system [^. These mechanisms might be 
independent with respect to each other or they might be mutually dependent. 
In case of Gompertz-type growth, bo = 0 and = 0 for n > 2. It indicates that 
the growth mechanism corresponding to the term bia is not dependent on the 
other growth mechanisms found in nature. In case of West-type allometry based 
biological growth process, bo — 0 and = 0 for n > 3. The growth mechanism 
corresponding to the terms bia and are simultaneously found in this type 

of system. Thus, the sole existence of the growth mechanism corresponding to 
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bia is reported and co-existence of growth mechanisms corresponding to the 
terms bia and 620 ^ are also found in nature. The sole existence of the growth 
mechanism corresponding to the term 620 ^ or a coexistence with the terms other 
than hia is still not reported or not even considered from the theoretical point 
of view in phenomenological approach. Here, such a possibility is explored from 
phenomenological point of view. 

In the phenomenological class represented hy (p = bo + 620 ^ with 5i = 0 and all 
= 0 for n > 2 , the analytical solution of y is expressed as, 

y = A'^ I cos[aT — 9] \'^ ( 8 ) 


Where, A = 0 = arctan/3, a = \/b^, P = Here, a may be treated 

as frequency of oscillation. 

The solution can also be achieved in terms of sine function as. 


y = 


A'^ I sm[ar-|-^] 


(9) 


Where, H = = cot ^ p. 

Though the solution can be expressed in terms of sine function as well as cosine 
function, while the cosine function would be considered to describe the be¬ 
haviour of the system in this communication. It is verified that the same fea¬ 
tures of the system can be extracted in both cases. 

The specific growth rate of this phenomenological class is expressed as. 


a = 



z) 


( 10 ) 


Where, 2 = Vbob^r. 

It is not possible to describe the behaviour of the system with the help of a single 
differential equation. It is because of the fact that y is not differentiable when 

[ar — 9] = . Therefore, the formation of segmented or piece-wise 

differential equations are possible that represent the behaviour of the system. 
The differential equation given below, 


1 <Py 

y dr^ 


b2-l^dy^ 


( 11 ) 


is valid for the interval (2n -I-1)| < [ar — 9] < (2n -1-3)^, where n = 0, 2,4,. 

. For the interval 0 < [ar — 0] < f and (2m -I- l)f < [ar — 9] < (2m -b 3)-| 
(where, m = 1,3,5,.), the differential equation is given by. 


1 d'^y 
y dr'z 


b2 - 1 .dy 2 ^ 
2 


( 12 ) 


For the phenomenological class corresponding to (/? = 60 + + ^ 20 ^ with all 

= 0 for n > 2 and satisfying the condition b^ < 46o&2; i-e; zeros of the 
function p is complex; the system is governed by the following expression. 


y = Pexp{—at) \ cos(wt 


7) 1 ^ 


(13) 
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Where. P = .r = r. = ^b„b, - ?, 

In this case, a = may be termed as damping factor of the oscillating system, 
and w = \lb0b2 —may be termed as angular frequency of the oscillation. 


3 Discussions 

In the proposed class of phenomenological universalities defined hy ip = b^ + 
620^; a, being a function of 5 o and 62, and A'^ can be treated as the angular 
frequency and amplitude of oscillation respectively. It is found that the am¬ 
plitude of oscillation of the system is affected by the variation of both 69 and 
62- In both cases, amplitude decreases with the increase of 69 or &2- In case 
of change in frequency of oscillation, 69 and 62 play same role, but the change 
in frequency is more affected by the change in 69 than the change in 62. The 
frequency of oscillation increases with the increase in 6g and 62- The change in 
nature of oscillation with the change in 69 and &2 is represented in figure © 
and ©• The nature of oscillation can also be characterized by the sharpness 
of oscillation that does not depend on 69. It is totally controlled by 62 (as it is 
observed in graphical representation). The sharpness of the characteristic curve 
increases with the decrease in 62. The change in sharpness of oscillation with 
the change in 62 is represented in figure Q. No such behaviour is observed 
in case of 69.This proposed phenomenological class might be useful to describe 
phenomenologically sustained oscillations observed in different physical systems 

[SI nil 131 ISO]. 

The change in specific growth rate with the change in 69 and 62 are represented 
in figure © and figure ©. It is found that the frequency of oscillation of spe¬ 
cific growth rate increases with the increase in 69 and &2. As a consequence, 
the time period corresponding to the change in specific growth rate decreases 
with the increase in the values of phenomenological coefficients. The frequency 
of oscillation of y is identical with the frequency of oscillation of specific growth 
rate, as expected from equation © and (HOD. 

The phenomenological class corresponding to p = b^ + bia + &2a^ with the 
condition bf < 4&9&2 leads to damped oscillatory nature of a growing system 
when bi and 62 are greater than zero.. The amplitude of oscillation is given 
by Pexp{—(JT). Therefore, the exponential decay in amplitude is governed by 
the magnitude of 61 and &2, as it is expected from equation da. The decay in 
amplitude increases with increase in bi and with the decrease in 62, as shown 
in the figure © and figure ©. The decay does not depend on bo as it is not 
function of 69. Thus, it can be concluded that the damping of the system de¬ 
pends on 61 and 62- It is practically governed by bi because of the fact that 
an undamped system may be considered for either 61 = 0 or 62 —>■ 00 but that 
is not feasible in reality. The magnitude of P increases with the increase in 
bi- But, the same decreases with the increase in 69 and 62. The frequency of 
the damped oscillatory system represented by this phenomenological class may 
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be represented by w, as it is found in equation (TT^ . It increases with the in¬ 
crease in bo and 62 (shown in figure ([ 6 ]) and figure (O) but the same decreases 
with the increase in bi (shown in figure (IS])).It is interesting to note that the 
damped frequency (w) of oscillation is less than the undamped frequency (a) 
of oscillation. The difference between a and w depends on bi. The amplitude 
of oscillation shows exponential growth for 61 < 0 and &2 > 0. This proposed 
phenomenological class might be helpful to describe oscillatory behaviour of a 
system accompanied by damping (may be positive or negative). 

The characteristic equation of undamped oscillation corresponding to the phe¬ 
nomenological class (fi = bo + 620 ^ could be obtained from the characteristic 
solution of the phenomenological class represented by = &o + bia + 620 ^ by 
considering bi =0. In a similar fashion, the undamped frequency (a) is ob¬ 
tained by considering 61 = 0 in the expression of damped frequency (w) which 
is less than the undamped frequency. Therefore, the change in frequency due to 
damping is solely determined by the term bi , that also controls the decay in am¬ 
plitude due to damping. Thus, it can be concluded that the mechanism involved 
in the system represented by the term bia is independent with respect to other 
mechanisms and solely responsible for damping. In absence of the terms bo and 
bi, the phenomenological class is represented by = 620 ^ and corresponding 
growth function is given by. 


y = ( 62 T + 1)'’2 (14) 

The same can also be derived from equation (IT3)) . In this case, 62 = I leads 
to a linear growth of the system. A system showing always linear growth with 
respect to time is not normally observed in nature. 

When the values of bo, bi and &2 are at the threshold condition bf < 46o&2j the 
oscillatory nature along with damping is about to cease, as represented in the 
figure ([ 8 ]). It is also found that the behaviour of the system is similar to the 
critically damped oscillatory nature of a classical oscillator when the the val¬ 
ues of the phenomenological coefficients are just above the threshold condition 
bi < 4 & 0 & 2 - This nature is represented in the figure ([5]). The behaviours of the 
damped oscillation and undamped oscillation of the phenomenological classes 
are very similar to the behaviours observed in case of a classical oscillations. In 
case of a classical oscillator, the undamped oscillation is characterised by the 
mass of the system, the stiffness factor or force constant and the damping deter¬ 
mined by the interaction of the system with the environment. The withdrawal 
of that interaction causing damping generates an undamped oscillation which is 
also found in this phenomenological approach. The damping of a classical oscil¬ 
lator depends on the property of the surrounding medium. Mass of a classical 
oscillator contributes in determining the damping and frequency of oscillation. 
Therefore, as an analogy, it can be concluded that the term bi plays the same 
role as the interaction of the system with the environment in case of classical 
oscillator. The terms, bo and 62 , determine undamped frequency of oscillation 
that is governed by force constant and mass of a classical oscillator. 
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4 Conclusions 


In the present communication, the oscillating behaviour of a growing system is 
considered in the framework of phenomenological approach. The variation in os¬ 
cillating behaviour is analyzed with the variation of different phenomenological 
coefficients and is represented graphically. The variation in specific growth rate 
is also studied. The phenomenological class representing damped oscillatory na¬ 
ture is identified with the conditional relationship between the phenomenological 
coefficients. In this connection, an analogy between the behaviour of damped 
classical oscillation and the damped oscillation in phenomenological approach 
is considered. The term solely responsible for damping is also pointed out. An 
effort is given to identify the physical significance of each term involved in the 
phenomenological description from this analogy. In a nutshell, the proposed 
phenomenological description is able to address the oscillation of a system with 
constant amplitude. The phenomenological classes may be useful to describe 
the behaviour of oscillatory physical systems [271 EH ESI El]. It may be benefi¬ 
cial to address oscillatory growth of a tumor. These proposed phenomenological 
classes might be helpful for the experimentalists to consider their dataset from 
a totally different point of view and to extract a suitable model for describing 
the temporal evolution of the system. 
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Figure 1: (Colour online) Plots of oscillatory nature of the variable y with bi = 0 
and 62 = 2. From top to the bottom the values of the parameter bo are 0.02, 
0.82 and 1.62. 



Figure 2: (Colour online) Plots of oscillatory nature of the variable y with 61 = 0 
and bo = 0.2. From top to the bottom the values of the parameter 62 are 0.42, 
1.42 and 2.42. 
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Figure 3: (Colour online) Curves of specific growth rate with 6i = 0 and &2 = 
0.4. From left to right the values of the parameter bo are 2.0 and 1.0. os 11 





Figure 4: (Colour online) Curves of specific growth rate with bo = 0.02 and 
bi = 0. From left to right the values of the parameter &2 are 3.6, 2.4 and 1.2. 
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Figure 5: (Colour online) Plots of damped-oscillatory nature of the variable y 
with bo = 4.0 and 62 = 0.4. From top to the bottom the values of the parameter 
bi are 0.8, 1.2 and 1.6. 



Figure 6: (Colour online) Plots of damped-oscillatory nature of the variable y 
with bi = 0.9 and 62 = 0.4. From top to the bottom the values of the parameter 
bo are 2.0, 4.0, 6.0 and 8.0. 
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Figure 7: (Colour online) Plots of damped-oscillatory nature of the variable y 
with bo = 2.0 and bi = 0.9. From top to the bottom the values of the parameter 
62 are 0.3, 0.4, 0.5 and 0.6. 
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Figure 8: (Colour online) Plots of over-damped nature of the variable y with 
bo = 2.0 and 62 = 0.4. From top to the bottom the values of the parameter bi 
are 1.6 and 1.8.. The condition defined by bg = 2.0, bi = 1.8 and 62 = 0.4 is 
just above the threshold of the condition bf < 4bob2 and represents a behaviour 
similar to overdamped nature of a classical oscillator. 
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